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A recurrent problem in industries is that of locating • 
new facilities and the selection of materials handling equip- 
ments for moving mterials between the various locations. The i 
problem can occur at various levels of coraplexities. In the ■ 
past, most of the investigators have tried to solve the problem, i 
of location and materials handling equipment selection as tw , 

independent problems. However, it is a well known fact that j 

the selected materials handling equipment influences the faci-. 
lity location and hence the layout of the shop.v^TUn the present : 
vsorkjii^^erefore,^ an attempt is made to solve the problem of 
facility location and handling systan selection in an integrated 
manner , sXtwo alternate approaches are proposed. Eist approach 
uses ZERO-ONE procedures where as in second approach, the prob- 
lem is structured as multistage sequential decision problem and 
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is solved using dyaamic programming procedure. A meiliod to ^ 

overcome the dimensionality problem in the dynamic programming ir-y 

¥ 

procedure has also been suggested.'^ 

,, A number of problans have been solved using botti the 
approaches. It is found that compulation time required to solve 
a given problem by dynamic programming procedure is much more 
than -that required by the Zero -one programming pmcedure. How- 
everj the extra information obtained in a single computation run 
of dynamic program outweighs the additional time required. ] 

( Both procedures give optimal solution, nevertheless 
there are limitations to each of the proposed methodologies, 
large memory location requirenent in dynamic programming proce- 
dure has restricted its application to problems of 5 to 6 flOY/ 
paths and two types of material handling systems that can be 
shared by two or more flow paths. However much larger problems 
can he solved using zero-one programming procedure. A problem 
involving 3 new machines, 9 flow paths and 4 alternate materials 
handling systems has been solved: using 0 - 1. programing proce- 
dure in around 30 minutes of corapu ter time on IBTI 7044/1401 
system. ^ 


I 



CHAPTER I 


INTROKJCTIOH AHD HTERATOEB REVIEW 

In the overall spectrum of manufacturing management, plant 
layout and materials handling plays a vital role. The productivity 
and "the cost of production of a company are very strong functions of 
its Ipeatf^j the layout and materials handling methodology. It is 
interesting to point out that for every ton of finished goods that go 
out of a factory door, atleast fifty tons must he lifted or carried 
around the plant to produce it, upto QQffo of all indirect labou^^ cost 
in a plant and nearly one half of the tpinl cost of manufacturing can 
be attributed to materials handling. To gauge the importance further, 
atleast one lupee out of every ten rupees the industry spends on plant 
improvement today goes to the purchase and maintenance of materials 
handling equipment. The above figures are enough to reflect the impor- 
tance of layout and ffiaterials handling and emphasise lb e S^ct that a 
saving of a small fraction of materials handling cost here and there is 
a significant total cost reduction. 

Materials handling andp lant layout are so dLosely interwoven 
that it is difficult to distinguish between them. It is impossible to 
install an effective handling system without considering the plant 
layout and vice-verca. One of 1he widely faced and a bit neglected 
class of layout prolJlem is the Revising and Inproving of Existing 
Plant Layout, To maintain operating efficiency or^er a period of time, 
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a given plant layout must ordinarily be modified and improved in some 
way. Iven in a newly constructed plant, management may find it nece- 
ssary to make revisions in botti the process and layout because of rapid 
technological advances and unforeseen changes in demand. The dynamic 
nature of competitive industry therefore, compels progressive manage- 
ment to make layout changes for various reasons. Some of these reasons 
are as follows : 

1. A plant Isyout may be revised (for eoq^ansion in production 
capacity) to meet a relatively permanent incr^se in the demand 
for its products. Such expansion ngiy involve pmviding faci- 
lities for the production of new products. This may require 
the addition of new facilities to the plant or the conversion 
of existing unused facilities to the output of new related 
products, 

2. A plant layout may be revised to attain greater competitive 
efficiency by adopting various tec±ino logical improvements 
that effect the plant and its operation. Such advances in- 
clude improvements in product design. 

3. Plant layout revisims may be the result of work -simplifica- 
tion programs that improve production processes and layout 
arr ang ement s . 

The most frequently securing layout modification is due to 
expansion in production capacity vhich requires the addition of new 
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facilities to the plant. Depending upon the situataon there mi^t be 
a considerable material flow between the new and existing facilities 
and also among the new facilities^ ia-volving significantly large 
materials handling cost. By making a judicious allocation of new 
machines and proper selection of material handling equipment, this 
materials handling cost can be minimized. 

Inspite of the fact that optimum location of facilities would 
yield considerable amount of cost reduction, only a few methodologies 
have been developed and reported in the literature. A review of the 
literature available in this area is presented in the following 
section. 

1.1 IiITERATDRE EBittEW 

The literature on the subject is sparse and in the past, 
investigators have primarily been concerned with the trial and error 
solutions , 2>. The first analytic technique was reported by 
ANDRE S. BIDDSOHEDIER and JAMES M. MOORE tC3> in 1961. They considered: 
the ccntinuous lo cation problem o f lo eating a single new machine in an 
existing plant. The proposed graphical technique consists of drawing 
level curves (iso - Gost Curves) for the aitire plant and choositig 
the suitable location on the minimum cost curve. The technique, how- 
ever, was restricted to only one new machine. 

WIMMEET m and COMAY 45) also examined this problem but 
limited thanselves to discrete candidate locations. They assumed 
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that the new Eiac±Lin.es will have no interaction, in terns of materials 
flow, w’ith the existing machines. This situation is certainly possible, 
but rather seldom found in practical cases. 

JAMES MOOEE (63 contradicted the WIMilMlT - CONlAY view, and 
developed assignment model which considers that all new machines will 
have materials handling contact with Ihe existing machines. Discrete 
candidate areas for locating the new machines were still used and no 
flow of material among the new machines was considered. The proposed 
technique uses ’’Hungarian Meihod for Assignment Problem" for locating 
the new machines* 

It needs to be pointed that in all the above reported lite- 
rature except Oonway's and Wimmert’s work, no flow of material between 
the new machines has been considered. TEANOIS 0., 8^ has tried to 
extend MOOEE' S level curve technique for this case but his analysis 
holds good only for one dimensional cases, i.e. , all the existing and 
nev/ machines are in a strai^t line. Moreover, the analysis turns out 
to be quite complex for more than two mew machines. This method can also 
be extended to two dimoasicnai case vhen the flow of material is equal 
in all flow paths. However, this is quite an im^sracti cal situation. 

Another work reported with flow among the new machines is 
due to HOSER G , m&IH and JACK D. EOGERS (9). They have given ^ 
analytical algorithm for Moores level curve technique for single new 
machine by using minimal property of the median. This algorilhm 
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assures optiml solution wliich may be. unique or lie within a range, 
For the case of locating several new facilities, Tergin and Rogers 
have instituted a computer search process comprising the following 
three major steps s 

1. Each nevi? machine in turn is located at a temporary 
optimum point assuining other m + n - 1 (m = number of 

* existing rja chines, n = number of new machines) machines 

fixed at their current locations, 

2. After all n new facilities are located in this manner, 
the process is repeated. 

3. Readjustnaat process is continued until no further move- 
ment occurs during a complete round of adjustment eva- 
luations. 

At ®ch step in the above process, the computer subroutine 
for single machine was used to find ihe temporary optimum location. 

The proposed search procedure is computationally fast but does not 
assure 1he optimal allocation. / 

TOUIBOUGrH 5^0^ in his unpublished M,S. work has given a 
dynamic programming approach for joint facility location and equip- 
ment selection problem. He suggested that materials handling cost 
very muda depends upon the type of materials handling equipment used 
for particular flow path, haice the problem of locating new machines 
and selecting suitable materials handling equipment for eada flow 
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path, should be considered sinultaneously. In structuring this problea 
as a multistage sequential decision problem, each new machine and the 
corresponding flow paths are considered as the stages while the number 
of potential locations a-yailable for the new machines are treated as 
state yariables. Material handling systems and material-handling costs 
are considered as decision variables and return function respectively. 
This approach gives optimal solution for the situations having no 
bu dge to ry constraint for the purchase of EO^terials handling equipments. 
Then the investment constraint is introduced, the algorithm gives sub- 
optimal solution Which sometimes might be quite far from the optimal 
solution (10). 

1.2 ERESBra WORK 

Realizing the importance of facility location with reference 
to cost reduction programme and recognizing the fact that no good 
approaches are available in the literature to handle practical facility 
location and materials handling problem in an integrated manner, the 
present work has been undertaken,. 

The present investigation is in line with Willibough' s : ^ 
work and two methodDlogies are proposed for finding overall optimun 
solution of Joint Pa cil tty location - Handling System selection prob- 
lem with budgetory constraint. 


OHAPTER II 


PROBIiM IlimilllOH AND SORMaLi^IIOH 

2.1 PROBLEM DEEINITIOF 

As mentioned in the previous ciiapter, the oh j active of the 
present work is to find an analytical solution procedure for a Joint 
Pacility Location - Handling Equipment Selection problem. The problem 
can more esqjlicitly be stated as follows ; 

"Eor an existing plant, the management has taken a decision 
to add some new machines, may be to meet permanent increase in demand* 
There is flow of materials between the new and existing facilities 
and also among the new machines. Eor every machine pair having mate- 
rials flow, there may be a number of alternate materials handling 
systems that can be used to transfer the materials. The objective is 
to determine simultaneously ihe optimal locations for the new machines 
and the optimal materials handling equipments to be used among the 
■various machine pa ir s so as to minimize Hie total materials handling : 
cost subjected to the condition that the investment on Ihe new mate- 
rials handling equipments does rot exceed the available finances." 

2.2 ASSDMPTIOFS 

The problem as stated is quite complex in its o-verall 
spectrum aid an analytical solution for this generalized case seems 
to be impossible. However, by making some judicious assumptions the 
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problem e an be simplified. The assumptions are introduced without 
affecting too Euch the practical features of the problem. The main 
assumptions are : 

1, lETlRffllTI STIC Flow OP miERIAI 

For the preseat aaalysis inter -facility materials flow 
is idealized to be deteministicj ■vh.ich is not strictly true. This 
assumption is introduced on the grounds that if the planning data are 
accurate, anall imbalances due to variable flow will resolve themsel- 
ves in the Icng run. 

2, DISCRETE LOCATIONS 

Facility location problau can be handled using either dis- 
crete or continuous locations. If a grid was to be superimposed over 
an existing plant layout, the new machines which themselves are con- 
sidered as point locations, cjou Id be located only at the intersecticsa 
of the grid lines in the candidate area. A continuous location prob- 
lem would exist if the new machine could be located ariyvh ere with in. 
the candidate area, i, e. , between the grid lines as well as on them. 
Whan a new plant is to be built, the. only layout restriction is the 
building configuration. Therefore, in such cases continuous locatj-ons 
make more sense than discrete candidate areas. However, when new 
machines are to be fitted "in to" the existing plant layout, cons- 
traints due to building configuration end present set up result 
in only a few discrete potential locations. 
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3. OTHER ASaJMmOlirS 

A few ofher logical assumptions are j 

(a) Materials handling costs are directly proportional to the 
distance moved. 

(b) The materials handling cost per unit of distance are the 
same for a particular handling system for all machiiie--pairs 
irrespective of the direction of nateral flow. 

(c) Between a machine-pair only one type of handling system is 
used. 

(d) Wages of the materials handling equipment operators are not 
included in the present fo3wnulation. 

2,3 EORMI/iTION : ZMO - OlE PEOGEii.M 

Eor the problem defined in Section 2.1 and with aforemen- 
tioned assumptions, two alternate approaches are proposed. This 
section deals with the ZERO - ONE programming formulation where as 
in Section 2,4 a djnamic programming approach is pres ait ed, ; 

2,5,1 Definitions 

(a) Plow PATH i Any facility-pair having materials flow 
between thau is termed as a flow path. 

(b) lOG/i-TION OOMBINATION j Any f easible assignment of new 
facilities to iiie potential looaticais is termed as a 
location combination. 
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2, 3.2 notations 

Number of new iaachla.es 
n Humber of poteatial locations. 

M Humber of existiag machines. 

ICS Total number cf possible location combinations. 

P Total number of new flow paths. 

I Suffix for flow path. 

J Suffix for alternate materials handling system. 

K Suffix for location combination. 

P Set of all J's. 

PI P1CP comprising those J’s for which corresponding 
materials handling system can be shared by two or 
more flow paths e.g., handling trucks. 

P2 P2CP comprising all J's for which the correspondlig 
loan dling system can not be shared by more than one 
flowpathv, e.g.j conveyors. 

■ A (l, J, k) Eraction of Jth (j Pi) type of handling equipment 
required for Ith flow path with Kth location combi- 
nation. 

Z(l,J, K) AO-1 variable which is : 

1 If Jth handling system is used for Ith flow 
path with Kih location - combination. 

0 Otherwise, 

C (l, J, K) Material handling cost incurred if Jth handling 
system is used for Ith flow path with Kth 


location combination 
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E (l, J, K) Purdaase cost of Jth handling systen (j £ P2) j 

! 

for Ith. flow path wilii Eth location combination. 1 

PO (j) Purchase cost of Jth handling system (j € Pi). 

S (j) Unused capacity of Jth type of handling system 

in the present layout. 

G Maximum investmaat fbr the purchase of handling 

e(3A.ipments. i 

y (j, e) Eumber of Jth tjrpe (j € Pi) of handling equip- 
ment required with Eth location combination. j 

2,3.3 Pormulation i 

ibr a particular location combination E, the -various i- 

constraint equations are ; I 

■■ '■■■-I 

1, HMPIIEG STSIM SEIEGSIOE COESTR/ilET 

It is desired to use only one type of handling system for i 

a particular flow path. Ma-thematicaliy, for each flow path, this | 

constraint can be expressed as i 

5Z. X (l, J, E) = 1 I I = 1 , E ( 2 . 1 ) 

:j 6 p ■ ■ ] 

2 . lEEBS®®!! GOHSTEAIIf j 

■ ' ■■■ .; ■ ^ i ., ■" .' ' " ■ ■... ■ : ■ '' '■ i '.■I 

, ■ ' ' ■ ■■ ■■ . , ■ ■■ ■ . :■ '.-'"i ■ '■ " ■. • .| 

J-th (j ^ Pi) type of handling system is al-ways bought J 

in # 10 le numbers, though the requirement may be in fraction?, e.g. , 

if 3.6 fork trucks are needed, then 4- (higher integerniyalue) trucks have 

: : ■■ ' , : .'■■ -O' ' ■ . , : / ■■ ■ j 

to be purchased. This can be acertained by adding following cons- | 

traint for Jth type of handling system (j & Pi). [ 
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F 

^ A (l, J, K) . X (I 5 J, K) - Y (j, K) :<S (j) ; all J , J € P1 
^ 1=1 

( 2 . 2 ) ' 

TWiere Y (j, K) is an integer and gives the number of Jth 
type of handling equipment required for Kth location combiaation. 

The investment constraint then can be written as t 

^ E (I,J,IC) . X (l,J,K) + ^ Y (jjK) . PC (J) ,^ & 

1=1 J t P2 JeP1 

(2.5) 

3. OBJECTIlffi HJNOTION 

The optimization criterion is to minimize total materials 
handling cost. The obiective function con therefore, be vTritten as : 

Min Z (iC) = 0 (i,j,k) . X (i,j,k) ( 2 . 4 ) 

p 1=1 

Equations 3,1 to 3.4 fona a pure integer program in which 
X (IjJjK) are 0 - 1 variables where os Y (J,K) can take any integer 
value. Y (JjK), hov^ever, can be represented as a aum of 0 - 1 Yaria- 
bles. To illustrate, suppose a variable y(j,K) has a deteioinable 
Upper bound U. Than, vherever Y (j,k) appears an equivalent binary 
representation can be substituted. 

Y (J,K) = 1 (J,K) + 2 oOg (J,K) + ... + (J,K) + 

: ... + (j, K) y 

T/ihere = 0 or 1 and value of 1 is choosen such 1hat 

. 0 . 


( 2 . 5 ) 
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Combiniag equation 3»5 with, equations 3»1 ho 3.4 a pure 
0-1 program is obtained for the Eth location combinat ion. There will 
be in all a total of 'IC such programs and the oTerail objective 
therefore, is to find that value of K and corresponding Z (l, J, E) 
for vtiich 1he total materials handling cost is minimum, i.e., 

( 2 . 6 ) 


Using the dynamic programming strategy to solve the faci- 
lity location - handling system, selection problem, requires that it 
first be recognized as a multistage, sequential decision problem. 

This implies that two. types of variables, namely state -variables end 
■the decision variables together vi-feh stages an d performance criterion 
havo-tn be first identified, 

STAG-ES : Each new facility : as well as each flow path is consi- 
dered as a stage, A new facility will be called a 
major s-tage gad the corresponding flow pa1hs are 
termed as ininor Stages.' 

STATE VARIABLES ; Eor simplici-ty in structuring -the problem 
and subsequent formulation, two directional state 
vectors are introduced, i.e., state vector Z in 
the horizontal direction and Yj in 1he vertical 


Z (e) 


Min 
E = 1 , IC 

2.4 BMAMIC RROGRAMBTG EORMUIATIOU 


2 . 4 . 1 Gen er al 


dir ection 
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Vector Z = ( • x^) is a (S + l) dimensional 

vector. The element represmts the finances available for the 
purchase of handling equipments and i = 1, 3 is the fiactdon of 
ith handling system available, 

= (y-t 5 Jr) gives the location combination of the 

Old 

J facilities, 

COUTROI OR lECISION VAEIABIES ; Alternate materials handling 
systems are considered as the decision variables, 

PERKBJMCE GRITEEttOR : Materials handling cost is taken as the 
perfoimance criterion i/iiiich determines the effectiveness of 
a given decision. 

2,4.2 notations 

I Suffix for the major stage. 

J Suffix for the minor stage 

S number of alt eina.te mterials handling systems, 

P Set of all materials handling systems, 

PI P1CP set of all materials handling systems iflfciGh can 

be shared by tvo or more flow paths, 

P2 P2CP set of all handling systems which can not be 

shared by two or Bore f].ow paths, 

U (l, j) Set of all admissible handling ^stons for Jth minor 
stage of the Ith major stage, 
u Decision variable (u € U), 
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PC (u, Yj, I, J) 


PC (u) 

mo (u, Yjj I, j) 


mm. (u, Yj, I, j) 




P (X, Y^, I, j) 


MHC (u, X, Yj, I, 


K (I) 
Ax 


Purchase cost of uth handling systsa for Jth 
minor stage of thelth major stage (uC P2). 
Purchase cost of the uth handling system (u 4: Pi) 
Material handling cost for the Jth minor stage 
of ihe Ith major stage if the uth material hand- 
ling systan is used with location vector Y^* 
Praction of the uth handling equipment required 
for the Jth minor stage of the Ith major stage, 
with location vector 1^4 

State transfoimation function for X and Y direc- 
tional states respectively. 

State vector in horizontal direction. 

State vector in vertical direction ibr the Ith 
major stage. 

Optimum handling cost for all the flow paths 
of (l - l) facilities and J flow paths of Ith ; 
facility with resource vector X and location : 
vector Yj. 

j) Materials handling cost for the Jth minor stage 
of the Ith major stage if the uth handling 
system is used with state vectors X and Y^ in 
horizontal and ver tical directions respectively . 
Humber of minor stages in Ith major stage, 

Sta-fce increment vector. 
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2,4 »3 Ibrmulation 


Using the notations of the previous section, the recursive 

relation can be written as : 

Min / 

P ^|MHC (u,X,Y^,I,j) + P 


I, J J > 1 


Mn 


(2.7) 


^ |MHG (u, X, Yj, I, J) + P (X), \ 

Gy (Yj, j), 1 - 1, K (I - t)) j if J = 1 

Ihe iteration begins by specification of the minimum cost 
function at the first rinor stage of first major stage. 


P (X, Y^, "l) =u^ (u, X, Y^, 1, l)| • 

MHO (u, X, Y^, I, j) for equations 2.7 and 2.8 can be 


( 2 . 8 ) 


stated as follows j ■ 

If u 6' P'1 

MHC (u, X, Yy, I, J) = MG (u/ I, J) If MBR (u, Y^, I, j)^oc 

Or 

■PG: (u):':^' \'; 

^ ^ ^ ^ ^ ^ otherwise (2.9) 

u 6 P2 

MHO (u, X, Yj, I, J) = MHD (u, Y^, I, j) If PO (u, Y^, I, U) x 

(2.10) 

=' 00 ■ \ Otherwise : 

State transformation function G^ and Gy are defined as ; 

Gy (Yy, j) = (y^, y^j yj.-]) If J = 1 and 

Yf (y -i J y ^9 ... 


Yt if J > 1 


(2.11) 
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If u t -then 

<■ Sjj (2) = (x^5 - IfflE (ujT^jIjj).. Xg) if lIHE (u,Y^,I, j)^ x^ 

” (^0 ~ ^-]j •••) — MHR (u 5 Yj |X y J ^ • ‘^g) 

if IHE (u,Yj,I,J) 2;^^ and PC (u) x^ 

= 0 otherwise. 

If u -f. P2 then 

~ ^^0 ~ ^1’ ’** 

= 0 otherwise, 

'0' is a null vector. 

Using equations 2.7 to 2.13s ^ dynamic cocipu tat ional proce- 
dure can be Set up. 

Both the formulations discussed above have been computerized, 
Iheir computational aspects esre discussed in following chap tar, 

2,5 miA KE THE 

Principal (^ta required for the analysis are : 

1, Equipmaat requiremait. 

2, Materials handling cost. 

3, Permissible location Combinations. 

To obtain these, following data have to be compiled from 
the various sources (from the existing plant, handling equipment 
suppliers etc, ) . 


(2.12) 


if PC (u) ^ x^ 

(2.13) 
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2,5.1 Data Compilation 

1. LOCAOION MATRIX (l) 

This matrix gives 1he alternate locations for the new faci- 
lities, There is a possibility that due to some technological reasons 
a particular nachine can not he located at some locations. The element 
of this mati'ix are either 0 or 1. 

D (l, J ) = 0 If Ith new nachine can not be located a.t 

location J, 

_ 1 If Ith new machine can be located at location J 

Prom this matrix possible location combinations can be 

generated, 

2. DISTANCE MATRIX (d) 

This matrix can be developed by the layout drawing of the_ 
plant. The physical distance between any two layout points can be 
measuxedj depending upon the conditions imposed by the material hand- 
ling system used and the present layout. Distance matrix is developed 
for two basic types of handling shipment movement namely straight 
line movaaent case and rectangular movement case. 

If . R = Humber of po^ 

M = Humber of existing machines. 

■Then ■' ■ ■ ;.' 

D (l, J, K) = Distance between the location pair I-J (or 
existing machine if J > H) for Kth 'type of 


movment case 
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K = 1 Eectangular movement case. 1 

K = 2 Straight line movement case. 

3. MATERIALS HAEELIF& EQaiPMUM SEECIEI GAT IDES 

The specifications of the various materials handling equip- i 

ment can be sought from the suppliers. These specifications include 
the capacity, the purchase cost etc. The operating cost for the ; 

equipment has to be calculated, based on costs of fuel, routine mainte- 
nance, repairs, maintenance overheads and depreciation of the equipment, | 

4. MATERlAiS HMiaiEG TOIIJME (MHV). "V;] 

The basic method to define the materials handling volume ; 

between two facilities is from-to chart or matrix. This matrix ^ows i. 

the inter— facility materials flow in unit loads per unit time for a ^ 

particular handling system. i’ 

■ ■ ' V ■ ■ . 

MHV (l, J, K) = Materials flov/ between I - J machine pair ( 

i 

for Ith handling system. ; !, 

5. MATERIALS HAEIiDING ALTERNATIVES (MHA^ ? 

Depending upon the type of materials flow, this laatrix | 

giye8 "bile alt exnate iio^terdals handliag equipments that can he used I 

' ■■ . ■: . '■■■ :'v: . ■,/. ■'■'.a';'', ■■ , '■ , , 

for transporting materials between a machine pair:. This matrix is | 

coded in three numbers viz., 0, 1 or 2 . This implies that ; | 

MELA (l, J, K) = 2 If there is no flow of material between I 

machine pair I - J. I 

■■ -'''V ^ 'I 

■. . V.'": i'' 1 1'' i" ■■ i ."..,4 ■/ -'I 

^ 'tz'-' :""\':::'z :::^:iz^ 
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= 1 If equipmaat K can be used for bran sporting 
material between machine pair I - J. 

= 0 If equipment K can not be used for machine 
pair I - J. 

6. STMmRD ilMB PER MOI® 

Standard operation time can be calculated by studying the 
basic operations during a move e.g.j for a fork tmck, these opera- 
tions are ; 

a. At start : 

lower fork onpty 
Raise fork loaded 

b« Run 

c. At the end of run j 
Raise fork loaded 
lower fork empty. 

lime for these moves can either be calculated by time 
study or by using standard graphs usually supplied by the manufac- 
turers of the equipment. 

To this operation time, must be added contingency factors 
time v^aich include time lost if . some restriction by other loads or 
obstruction reduces the running speed. This time depends upon the 
operating conditions in the plant. Thus : 
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^ lormal Time 
per move 


^ Run time ^ 


lower and 
lift time 


Contingency ' 
factors time 


(2.14) 


To arrive at the standard time, allowances for personal 
time, fatigue and unavoidable delays are added to normal time. 


) Standard Time i 
I Per move j 

2.5.2 Data Analysis 


. ■) 

ime V 


formal f ime y 
Per mo Ye J 


Allowances J (2.15) 


Prom the data of the preYious section 5 e(j;.ipm oat requirement 
and materials handling cost can now be calculated. 


1. ESTimTIFG- EQUIMMl RECPIREi}/® 


Based upon the standard time and material handling ■volume 
from from~to charts a deterministic estimate for materials handling 
equipment requirement for a particular flow path can be made by using 


the following formula ; 


«HH (I, J. K) = 

TOiere AMHR (l, J, E) = Number of Kth type of equipment required 

along flow path I - J. 

MHY (l, J, K) = Material flow per monlii in unit load 

per unit tiae for Kth handling tan. 

ST (l, J, k) - Standard time in minutes, per move for 
Kth equipment between machine pair I-J . 
H = Standard hours per month. 

K = Utilization factor of Kth equipment. 
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2. MATERIAIS-HAlIBIiIHG COST 

Materials handling cost for transporting material "between 
a machine pair for a particular handling system is 

c (I, j, ic) = I) (I, j, k) * OG (k) *mY (I, j, k) (2.17) 

mere C (l, J, k) = Materials handling cost for transporting 

material between machine pair I - J with 
Zth handling system. 

OG (k) = Operating cost for Kth type of handling 

system. 



CHAPTER III 


COMPUTATIOmi ASPECTS 

In the pre'vious chapter two alternate formulations were 
proposed, solving them by conventional methods on computer mil require 
lot of computation tine and unmanageable high speed maaoiry locations* 

In this chapter computational aspects of the two approaches are dis- 
cussed. Section 3.1 deals with modifications in zero -one programming 
procedure to save computation time vliere as in Section 3.2 a method 
to overcome dimensionality problem in the dynamic programming approach 
has been suggested. 

3.1 MODIEICATIOI II ZERO-OIE program 

3.1.1 Solving Procedure Eor General 
Zero-One Poxigram 

At presaat several methods are available for solving linear 
programs of the 0 - 1 type. Best kiown among them are R ,B. GrOMOEY* S £ 1 1] 
algorithm for aolvtng linear programs with integer variables. He uses 
the dual simplex method and imposes the integer conditions by adding 
automatically generated new constraints to the original constraint 
set. 

Another algorithm for integer linear programs, developed 
by A.H. laid and A.G. DOIG <1123 , starts with a non-integer optimal 
solution and then finds the optimal integer solution through parallel 
shift of the objective function - hyper plane. 
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All the above algorithms require solution of corresponding 
ordinary linear program (without the zero -one constraint) for each 
iteration and therefore, are affected by round off errors. Moreover, 
these algorithms do not take fill advantage of the 0 - 1 property of 
the variables. lY/JJESCU 115} has suggested a different approach based 
on Boolean algebra concepts. Another approach has been developed by 
BAIAS 114J .. Both diiese methods attack directly the 0-1 program 
and there is no round off errors as operations consist of only addi- 
tions and subtractions. 

Many more procedures have been developed and reported in 
literature to handle particular type of 0 - 1 programs. However, for 
the present work Bala s' algorithm is used because of its generalized 
nature and it works quite efficiently for the type of problems in 
which most of the variables take 'O' values 1143 . Bor liiis study 
package program for modified Bales algorithm givac by MIZE and 
KHESTES y53 is used.^ 4 ; 

3 , 1 , 2 : Mo dif i ca t ion s 1*0 r Pr e sent Work 

The problaa in its generalized form can be stated as ; 

(3.1) 

with an overall objective function 

(z^) (3.2) 

Viftiere IC = Humber of location combinations. 




■»c 
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The most ob-vious procedure is to solve each program using 
Bales algorithm and select the combination with lowest -value of 
But computationally this is not economical because of the large 
computation time required to SDlve a 0 - 1 program. 

Another procedure is to couple all the (ic) programs to 
a single 0-1 program. However, not much saving in the computation 
time was achieved and even in some cases the time required was more 
than the total time required for solving the problems individually. 


Coming back: to tae generalized case, the problem can be 

physically presented for -the two variable case as win in Figure 

1, 2, 3 are -the convex feasible regions 
correspohding to different location 
combinations. .Exploration of each 
feasible region involves "the solving 
of a separate problem. /When problem 1 
is solved, a minimum value, 2 ^^^ , of* 
the objective function is acjhieved. 
Before solving problem 2, an additional 
constraint 

C„ X : ^ : Z / (3»3) 

^ 2 min . 

is intro diced. Depending upon the 
vector 0^, the constraint can take either direction A or B or 0. 

Cases A and B show that modified problem - 2 has a feasible, solution 
without 0 — 1 constraint. In case A, 0 — 1 feasible solution is 



HCUHE -> 3.1 GEAPHICAl 
REPEESElITATIOl]' OF ZERO-ONE 
EROGEAMMnra KEMJIATIOH. 
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present whereas in Case B there is no feasible 0 - 1 solution. In 
case 0, there is no feasible solution. Lot of computation time can 
be saved if it can be found before handj whether the present location 
combination is going to improve the obiective value or not. Case 0 
can be detected by coupling a revised simplex program with Bales 
algorithm package and using phase I to find the feasibility of the 
problem. The general solution procedure can nov; be summarized in 
the following steps. 

1» Bet (x) = Minimum value of objective function before I th 

location combination is tested. 

2. Add the constraint Cj Z ^ (1) in the Ith problem set and 

fiiad using phase I of revised simplex procedure, whether feasible 
solution exists or not, 

3. If no feasible solution is present, then, C I -f l) = Z^ ^ (l) . 

Take the next problem set and repeat the procedure. 

4. If feasible solution exists, use Balas algorithm for finding .0 - 1 
integer solution of the problem giving -value of ob jective function 
as Zj. Zj will be a high number if there is no feasible 0-1 
solution. 

5. Put ( I + l) = Min I Zj, Z ^^ (l)j . Repeat the procedure 

with next location combination. 

A logic flow diagram for the above procedure is given in 


Appendix - 1. 
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3.2 DIlffiKSIOM-iililTY EROBIEM IH 
HIOGMMMIITG- APEROACH 

Solviag Idle dynamic programming formulation of the location 
prohlem using conventional dynamic pxogramming computational procedure 
proposed by EIGffiiED BEMiU |163 is impossible because of the enormous 
computation time required and due to what Bellman Q\S)t calls ''the 
curse of dimensionality" - i.e., the excessively high ^eed memory 
locations required. Eor a problem consisting of two handling-systaas, 

5 location combinations and a state incranent of .1, the memoiy require- 
ment is of the order of 200 K which is beyond the capacity of mm^^ 
presoxt day commercially available digital computers. ROBERT E. lARSOH 
has proposed a state increment dynamic procedure to overcome the 
dimensionality problem. The technique consists of doing computation xn 
small bloclcs, but it is assumed that the next state after a control 
action is tahen, always lies in the neighbourhood of pres state. 

Por the present work this assumption is not valid. Therefor Sj a 

ble state procedure is proposed.* This new procedure has computati^^^ 

requiremmts that are signifimntly le the conventional 

case. requirement is alv^ys redu- 

ced, generally by orders of magnitude. 

^^This procedure is named variable state ctmaaic procedure, because 
the effective states in this procedure are blocks instead of grid 
' points and these blocks vary at each stage. 
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5.2.1 V/OlI/iBLE SW.TB DYmilC mOGR/JMMICJ& 

Philosophy behind the vex iable state dynamic procedure is 
that fbr a problem in which there are few stages and discrete decision 
variables, many neighbouring states have the same optimal solution. 

So instead of storing control action and objective ■value at each grid 
point o'f a multidimensional state field, only "the blocks ha-ving the 
same optimal solution and corresponding control action and objective 
value can be stored. The memory requirement can -therefore be reduced 
drastically, Oompu tat ion time can also be reduced by using ■the fact 
that if "the diagonal points of a block, one nearest and another far- 
thest to the origin (j 2 f state), have the same optimal solution, then 
every grid point with in -the block has the same optimal solution as 
the diagonal points. Hence, by using proper search procedure, so In tr on 
at many grid points can be found implicitly without actual computa-tion, 

Por the present work, convoational dynamic programing 

procedure is used for states ifi Y dir e ctfcn, vdiereas -variable s-tate 
dynamic procedure is applied for states in X direction. The variable 
state procedure is summarized in -the following steps : 

STEP s 1 

Discretize the states Z (j) wi-th incraaent A ^ 

.plained to section 2.3, 2 (j) iB = (S h- ,) aimaaaio»l vect=r and 
Hencs the dieoietiaed vames of X will fom a (S t 1 ) dtaensdcnal 
grid. Tteae dlsoretioed taluee of % are the points *ere oomputatlons 
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are to be made. The grid points are Identified by using: coordinate 
system with state (o, 0, ... o) as the origin. State! (j) cnn be ; i 

represented as Z (j) = coordinates of i 

X (j) are ( p j) , 

STEP ; 2 : 

Create a master list whi«±L stores the coordinates of 
diagonal points of candidate blocks. Initially at each stage this 
will contain a block having all the grid points within it. A block 
is identified by its two diagonal points, one nearest to and another E 

farthest from the origin. The coordinates of the two diagonal points [ 

are a (j), J =1, 2, (S + l) and C (j) , J = 1, 2, .. . S + 1 
respectively, | 

STEP : 5 I 

. ■ ■ ■ ^ ■ f 

If master list is mpiy - go to step 13. I 

■■■ ■ ' ■' '. . ^ ' ' .• i 

STEP S 4 ■ I 

. E ^ ^ . ■ ' ' . ■ I 

Take out first candidate block from the master list and ; 

put X (j) = a (j) #iere 3 = 1, 2, ... S + 1. Find Op tinal solution j 

ISO! at this point, using recursive relations of section 2,3 | 

STEP : 5 ".| 

Put X (j) = X (J) t. n h Xjl forall J. Wiere n = 2^ 

(k is any integer). ; 

If all X (j) < P (J)j, go to step 7 otherwise go to step 8, s 


STEP : 6 
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STEP ; 7 

Piad optinsl SDlutioa at point Z (j). If solution, is same 
as ISOIij go to step 5. Ottierwise go to step 8, 

STEP *. 8 


If n equal to 1, go to step 9* Otherwise put n = n/2 
and start ttie search procedure from step 5 with the previous grid 
point, 

STEP ! 9 

Incrffiient, individual X (j) Vs, as hSir from the origin as 
possible, such that optimum solution is the same as 1301. 

Diagonal points a (j)'s and X (J)' s give a block, having 
the same optimal solution ISOD. 

STEP .1 10 ' 

Hemove the block from the main block and divide remaining 
Space In pxoper candidate blocks, A partitioning procedure for 
different position of solution block in three-dimaasional grid is 
shown in Pig. 3,2, 


STEP s 11 


Enter the partitioned blocks in master list. 


STEP ; 12 


STEP t 13 


Search procedure is cmpleted, proceed to next dynamic 


programming stage* 
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iJie solution procedure is further illustrated ia figures 
3»3 io 3«5 hy solving a aiiall location problen with tv?o stages (flow 
paths) and one location cottibination. A case is considered where a 
new machine is to be added to a plant. This would introduce two new 
flow paths into the system. Sbr each flow path, either a handling 
truok (Decision - l) or a conveyor can be used for materials handling. 
Using dynamic piogrammiag procedure, handling equipment for each flovr 
path has to bo selected. Data for the problem are : 

Flow Path - 1 

Decision Yariable 1 2 

Equipment Requirement ,7 1; / h 

Handling Go 8 1 Rs, 500 : Rs, 400^^ h 

Purchase Cost of 

handling equipment Rs. 13,000 Rs. 19,000 

. Flow'Path - 2: 

Decision Iferiable . 1 : 2 / 

Equipment Requirement .6 1 

Handling Cost Rs. 595 Rst 500 

Purchase cost of 

hmdllng equipment R& 13,000 Ra 22,000 

Maximum finance available = 30,000. 

In ■this case S = 1, hence discretized states form a two 
dimensional grid as ^own in Figure 3 . 3 . Increment in direction - 1 
(representing finance available) is R$. 1,000 where as along dixecti.on-2 
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(representing fraction of truck available), it is . 1 . Using the search, 
procedure explained earlier and the recursive relations of section 2.3 
the state field is divided into number of blocks having the same 
solution. Pi^jure 3.4 illustrates search procedure for stage - 1 vhere 
as stage - 2 is explain ei in Figure 3 . 5 . lumbers within the block 
are the optimal control action for the two flow paths. The first 
number corresponds to the first How path vhile the second number 
corresponds to the second flow path. The symbols - 0 r and ft represent 
grid points where computation is performed, and the diagonal points 
of the blocks having the same solution respectively. For this example 
the optimal solution occurs at point . (30^000,0) , i.e. , use truck for 
both flowpaths. 

It is interesting to note taat if conventional dynamic 
pTOgremmtng would have been used, then control actions and objective 
values at a total of 600 states will have to be stored vh ere as using 
variable state dynamic procedure this number is as low as 22 . More- 
over computation is performed at only around 100 grid points as com- 
pared to 600 in conventional dynamic programming procedure, 

A computer program for this procedure has been written on 
FORTRAN - If. A general logic flow diagram for the program is given 
in Appendix - 2 , 


CHAPTER IV 


RESanDSAMD PISCUSSIOK 

In all, eight randomly gai era ted problems of various sizes 
are solved using both the procedures. Por problem - 1 data concerning 
acceptable locations of new machines, distances between tvro machines, 
material flow between a madiine pair and tie alternate handling equip- 
ment used, are borrowed from TOIIiIBOUGH' S paper £.10^; handling equip- 
ment specifications are taken from' catalogues obtained from M/s MdJally 
Bird Ltd and the standard times per move of various handling systems 
between a machine - pair are generated using random number tables, 

Eor this problem sensitivity analysis is performed for each parameter, , 
Other problems of varied sizes are solved to estimate the. computational 
efficiency and limitations of the two procedures, 

4.1 LEST ERCBLEM - 1- 

4.1.1 Problem Statement 

It is desired to locate three new facilities on a plant 

floor. Pour; existing faciliti es are 
already in operation in tie plant and 
their materials handlmg loading and , 
unloa ding stations are designat ed as 
Pj., Pjj, Pjjj and Pj^ as shown in 
Pigure 4,1, 

There are four locations on 

PIGURE 4,1 5 LCXiATIOIS OP EXISC- plant floor which may accommodate 
ING MACHINES ATO AVAILABLE) 

AEXAS, 
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on© or mors of fh.© facilities. liiese locations are designated "by 
Arabic Numerals in fi^re - 4,1« The other data are as follows : 

4.1.2 Input Data 

(a) Acceptable locations ; 

Facilities to be located are limited to the available loca 
tions as shown in Table 4,1, 

TABLE 4.1 

ACCEPTABLE LOCATIONS FOR TBE NEf FACILITIES - PROBLEM - 1 


NEW L_ 


LOCATION MJMBER 


MACHINE 1 

1 

2 

3 

4 

A 

1 

"'tv'.- 

V -I : 

^ ^ 1 

B 

1 

0 

' 1 - 

0 • 

0 

0 

1 

; 1 : 

h ■ 0 " ■ 


'1' at (a, Li) denotes that A can be located at LI 
where as 'o' at (B, L2) shows that B can not be located at L2, 

(b) Distance Matrix 

Equivalent distance between existing facilities and new 
locations for type of movanent followed by handling equipment is 
givoa in Table 4.2. 
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TiffilE 4.2 

msm^GisBwmm iwATiass - probim - i 


lOOATIOH PAIR 1 

mSMOE 

1 B METERS 

' MOVMEm CASE 

I STEAISHT' IBrE ' 
f MOVEMMT CASE 

1-2 

20 

: 15 : 

- 5 

18 

10 

- 4 

32 

23 ; 

2-3 

18 

13 V' 

- 4 

13 

12 

3-4 

32 

. 20 

PI - 1 

17 

: ^13 

- 2' 

20 

:.12' 

- 3 7 

: . 25 

17 

-4 

, 32 : 

■ 22 , 

PII - 1 

' ;'32; , 

20.:'' 

'"-'2' ■■ 

13.:..::. 

: 8 

-^3' ■ 

32 ■ : ^ 

■ ZO- 

- 4 

18 

IC 

pni - 1 

35 

25 

■ ; - P'' ; 

25 

17 

• 3 ■ 

22 

15 

- 4 

20 

15 

prv - 1 


25 



a:':; 


22 
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(c) Handling Equipment Specifications 

_ Specif icaticaas of alternate material handling systems are 

summarized in lahle 4.5. 

liBIE 4.3 

SPECIE! Oil lOHS OE MilERIAlS HAmLWG STSTMS PEOBIEM - 1 


SYSTEM I system! 0/^ACITY I pDROHASE JoPEEiaiNG 1 AVG. 

I CODE ! I COST X COST I SPEEP 


Hi-Stadc* 








(Electric 
Eork-lif t 
Truck) 

P 

2 Tons 

Rs. 

60,000 

Rs. 

1.05/ISit 

6 Em/h.r 

Mazdur-1015* 

(Hand-l‘ruck) 

Q 

1 Ton 

Rs. 

5,000 

RSi, 

,60/Mt 

2 Km/hr 

Bolt Conveyor 

s 


Rs. 

800/Mt. 

Rsi 

.45/Ton/^t 


Trolley 

Conveyor 

T 


Rs. 

1,200/iat. 

Rs. 

.55/Ton/l!4t 



*Trade names of Tincks manuifiictured hy MdKFally BSLrd & Go. ; 

(d) Standard Time Per Move s 

In Table 4.4 for each location pair the standard time 
per mo've for two types of tiucdcs under consideration are given. 
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TABIE 4.4 

STMIMD TIME PEH MOVE - PROEM - f 
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(e) Material Handling Alternatives 

The alternate handling systems vtoich can be used for parti- 
cular flow path are listed in Table 4.5. (p, Q, S) at (A, b) shows 

that for flow path A - B any of handling equipment Pj Q or S can be 
used for transporting material where as 'Z' at A - G represents.no 
flow of material between the machine pair. 

TABLE 4.5 

ALTIRNATE MiTEEI/il Hf.EI3LIIG SrSTEMS - PROBBM - 1 



A 

B 

0 

n 

PII 

illl 

EIV 

A 

♦ Ip'-* 

P,Q,S 

2 

2 

P,Q 

2 

P,Q 

B 

» f « 


Q,T 

1P,Q 

P,Q^T 

P,T 

: 2 

0 

•'■t * ■ 

* 4 .$: 

, f • * 

P,S 

: 2 

' 2 ; : 

P,QsS 


(f) Monthly Unit load 

Table 4.6 gives the monliily material flow in unit load of 
handling ©(jiipment for a machine pair. 

TABIE - 4.6 

MOITHLY mTERIAL EDOW IN UNIT LO/^BS PER MONTH 


HMEblNO j 
EQUIBCERT! 




PLOW 

PAQH 





A-B 

"J-EEI ' 

' A-HV 

giG ■ : 

B-FI 

B-fll 

B-pni 

C-EL 

C-EIY 

P 

100 

200 

125 

X 

125 

75 

225 

50 

25 

Q. 

200 

400 

250 

300 

250 

150 

X 

X 

50 

S 

200 




mXBM 



100 

50 

T 

r:'" X 


i'-iiSliii 

300 


'150 

450 

X 

X 
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4*#*1#3 B.0SuX"b 

The above mentioned problem is solved by using 0 — 1 pro- 
gramming procedure and final results are tabulated in Table 4.7. The 
value of the equipment utilization factor is assumed to be 0.7. The 
computation time required was of the order of 30 Min, on IBM 7044/1401 
computer. 


TABIE 4.7 

OlBBJiUL OPTIMQM PiiOIIiITI IOOATION - HAEIIIDrG SYSTEM 

siaiEaTioir policy - erqbibm - i 


PACILITY 

i Assicsaffii) ]i 

1 LOCATIOI 1 

BLOW 

PATH 

”1 SYSTEM 1 
iSELBGTED 1 

MONTHLY 
HfiNDLING 
COST HRs* 

T 

X 

CAPITAL 

COST 

Rs. 

A 

L2 

A-PII 

Q 

3,120 


5,000 



A-HY 

Q 

3,450 


5,000 

B 

L1 


S 

1,350 


12,000 



B-PI 

Q 

2,550 





B-PII 

T 

1,640 


24,000 



B-HII 

T 

6,200 


30,000 

0 

L3 

O-B 

T 

1,650 


12,000 



0-PIY 

S 

225 


8,000 



C-H 


765 


13,600 


Total Rs. 20.950 Rs. 1,09,600 
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In fable 4,7, corresponding io the flow path B-Et no 
anvestmont has been indicated. This is due to the fact mat it can 
share the trucks purdiased for flow path A-BII and A-HY. 

4.1,4 Sensitivity Analysis 

Sensitivity analysis has been performed for "the following 
parameters s 

1, Handling equipment utilization factor. 

2, Cost parameters. 

3, Constraint on Investment, 

(a) Handling Equipment Utilization Eactor 

Ibr determining the equipment requirements for each flow path 
a utilization factor K is introduced to take into account the break- 
downs and other unavoidable, delays. Geneiully ; for design purposes the 
•value of K is taken equal to ,7. However, an analysis is performed 
here for various values of K to detemine the response of optimal , 

solution to changes in utilization factor of the handling equipment, 
fable 4,8 gives the solution for various values of K. 
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T ABIE 4.8 

OVEEAltl OHIIMUM EACIIIiy ECXl/iTlON - HAlilllLIH'G- SYSTEM SEIEOTIOF 
POHCY WITH K AS A PARAlilETER - mQBIEM - 1 


ASSIGNED 

LOCATION 


K 

A 

B " 

"”o 

»6 

12 

11 

13 

.7 

12 

LI 

13 

.8 

12 

LI 

L3 


SYSTEM SELECTED 

A~EII A-EIV B-A B“II B~EII B-EIH 0-B C-flV 0-EI 
« Q S Q T T I S S 

Q Q S Q T I T S S 

Q Q S Q T I T S S 


MATEBIAL 
HANDUNG GOST 
Rs. 

20,950 

20,950 

20,950 


Results show that with the present set of data even if -the 
performance of handling equipments deteriorates to .6 or improves to .8, 
the optimum solution is not affected. 

(b) Cost Parameters 

r , 

Handling cost is quite nebulous in character because it involves 
human judgement and error. Therefore, a sensitivity analysis on handling 
cost to study Its influance on the various decisions is very desirable. 

The optimal solution obtained on varying each cost parameter indepoa- 
dently by 10^ m either side are given in Table - 4.9. 
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TABLE 4.9 


OVmALIi OPTIMUM mCIEITY LCXIlTIOI - MNDIIW SYOTM SELB3TI® 
POLICY WITH OPBR/VTING COSTS AS PAEiiMETER - PROBLEM - 1 


Normal 

cost 

Parametecs 


Vari at ion la Cost Parame ter s 



' A L2 L2 L2 L2 L2 L2 L2 L2 L2 

ASSIGNEE B L1 L1 LI LI L1 LI 11 L1 LI 

LOCmON 

0 L3 13 L3 L3 L3 L3 L3 B3 L3 


1 A-EI Q Q Q 

A~EIV Q Q Q 

B-A S S S 

B-Kt Q Q Q 

PLOW : B»HI T T T 

PATH 

B-pni T I 2! 

; O-B T : T T 

: C-PIV S S S 

^ •O'-EI. ■: V s. . : B ..S:,:; 

MATlBIi.L 
HANDLING COST 

IN Rs. 20950 20950 20950 


Q Q Q Q Q Q I 

Q Q Q Q: Q Q : 

s s s s s s ^ 

Q Q Q : Q. Q ■ Q f 

: ' ■ ; ■ ■ ■ ■ I 

T I T T T T I 

rp I 1 T T T j 

■ " ■■■ ■ ' ' I 

m I I T T T ^ 

'***.' I 

S s s s S:;;: S;,^ 

s s s s s s I 

■' f. '- '■■■ / ■ ■-:■■■ "" ''^l 

21663 20037 21185 20715 21899 20001 { 


Heaxlta giv® to ?!able 4,9 show that the optimal elution does 
not change wito toe 10^ Y^l&t#?h 4^^ cost component i.e., even if opera- 
ttog costs are calculated wronglyt the optimal solution is »ot influ 

oed^ aiieast' 
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(c) Investsaeat Gonstraiat 

Althoqgli tovestment is a management decision, yet it is 
desirable to study tiie effect of tlais parameter on o'verall optimal 
solution. If there is a significant saving by little extra in-vestment, 

"this altemati-ve can always be recommended for iranagement' s considera- 
tion, Por the present problem optimal solution is obtained for -various 
leyels of investment. The results are summarized in Table 4.10. 

TJffiLE 4.10 

OTBRAIIi OPTIMJM PACIIITY lOOiiTIOI - miDHIfG SYSTM SELEGTIOl 
MIOY TOH in/BSTMMT AS PARAMBTIR - EEOBIEIil - 1 

Investment Assigned System Selected Material 

Rs. Locations Handling 

A B 0 A-HI A-PIV B-A B-EI B-PIl B-Pni C-B G-PIY C-H Cost 

1.10.000 12 11 13 Q Q S Q . 

1.20.000 12 11 13 Q Q S Q T T T S S 20,950 

1.40.000 12 11 13 P P S P T T T P P 20,694 

1.50.000 12 11 13 P P S P T T T S P 20,347 

1.60.000 12 11 13 P PSPT T TS S 19,800 

Table 4.10 shows -that with the increase in Investment level, 
monthly materials handling cost reduces. It is for managemmt to decide 
whe-ther 1h.e saving is worth -the extra inves-toait or not. 
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4*2 ZmO-K)HE PROGmiMIlIG Vs. OTiMIO 
mOQBMmL'SGr PROCBIUEI] 

A'btieaiip'fc vvas also made iio solve Problem. — 1 by dynamic 
programming procedure but it failed because bf high memory require- 
ment, Therefore, to assess the applicability and efficiency of 
this procedure, problems of varying size were solved. Details of 
these problems are given in Appendix - 3* Appendix - 4 gives a 
sample output of dynamic programming procedure. 


The respective time required by both procedures for 
various problans are listed in Table 4. 11. As expected, time re- 
quired by dynaimto programming procedure is significantly more than 
that required by Zero -one programming approach. However, much 
more information obtained by dynamic procedure outweighs the ex- 
tra time reqiilred. In a single computation run, dynamic progra- 


mming prooedire gives optimum handling system selection policy for 


each location combination and for different levels of investment 
viaere as zero-one programming procedure gives only the .overall 
optimum policy for a particular investment level. 

The. "Obtained in ;’d3fcMie programmi^ 

procedure- ■■avoid ttoa.- extra \ccm®utaticp.- run 

n®@l©a ' 

tabla' iadhs" :fbr" 'lihfe' /Ca^es 
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If me tries to obtaia all the mformation vshich gets 
autoraaticBlly generated by the dynamic programming procedaro, 
using zero-one prograramlng , the computation time required would 
have been much more than the corresponding time for dynamic pro- 
gramming proceduro* 


TABLE 4.11 

zmo - OHE EROGEAMMING BROCEIUEE Y/S DICfAlIC 
EEOGlUaraG APBROAOH 


Compu ta tio n T ime on 

Problem Ho. of No. of Ho. of IBM 7044/1401 


No. 

New Machines 

Locations 

Plow Pa1hs 

0-1 Programming Dynamic Pro- 
Procedure gramming 

Procedure 

1 

3 . 

: 4 


3fMin. 

Oan not 
solved because 
of large memory 
locations re qd'. 

2 

■■ 1 ■ 

: 2 

2 

6 Sec* 

: 10 Sec* ’ 

3. 

■■ 2 ■ 

: ; 3' 

3 

8 Sec, 

: 14; Sec* ^ ' 

4 

: 2 


4 

1 Min. 11 

Sec. 2 Min, 40 Sec. 

■5' 

^ ■ 2; 

ii 

5 

1 Min. 1 

Sec. 2 Min, 10 See. 


:2 : 


i"' 

® f! 

. 2 Min. 41 

Sec. 6 Min. 46 Sec* 

7 

'"’2 - 

'‘S' 


■ 2 Min. 15 

Sec. 7 Min. 28 Sec. 

8 

:3' 

’•’ir 

£ 

0,-: 


Can not bo 
solved because 
of large memory- 
locations roqd‘. 
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4,3 LIMITATIONS OF PEOPOSEL API'ROACHES 


Both zero-one progranming and dirnamic progi«ming 
procedures give optimal solution to joint Facility location - 
Handling System selection pieblem, Nevertheless, there are limi- 
tations to those models. The models are limited to deterministic 
flow data* If flow data can be taken with a high degree of con- 
fidonce, then the prooodures should work well. However, in prac- 
tical situations it is more like iy that flow data will be stochas- 
tic and subject to random variations. Further stady is needed to 
determine the sensitivity of the layout solution to flow variations 

Another limitation occurs due to the size of problem 

that 3&an b© handled. Large memory requirements have restricted the 
appiioation of dynamio programming procedure to only upto 6 flow 
paldns, RcoblOTS up to 3 new machines, 9 flow paths have been sol- 
ved using 0 - 1 programming approach in less than 30 Minutes on 
IBM 7044/1401 systca, but computation time required increases ex- 
ponentially with the inpreaso in the size of problem. 


4.4 OONOLUSIC® 

V Iai^i^^iidj^3^;!e^ahQve' .mraatiohsd limitslfichs 

jTOmls© to ■ iln:v:«§|y :slaO 

problems require addition of 


lane,; 


5 to 4 tim 






titro duces less than 10 new flow paths, 
blems up to 6 flow pat 
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or z(|ro-one progromtog procedure (new flow paths nore than 6) 
can very well he used. However for larger prohleas, still 
better approach noodsto be developed. 

4.4 SUGGESa?IONS 3?0H HTRTH2E WORK 

-ft. few suggestions for the extension of the present work 
arc listed, below s 

1. In all the above reported literature, including present work, 
the objective was to nininize the materials handling cost and 
no woightage was given to ihe corresponding investnent. There- 
for© a better objective function could have been one, which 
simltanoously considers both investment as well as materials 
handling cost* One of 1he approaches may be to minimize the 
present worth of investment and materials handling cost over a 
predefined design period . Ihis will slightly nodity the o b- 
jectlv© function in case of zero-one programming procedure but 
to take into account this aspect in dynamic programming formu- 
lation can bo a dhallmglng pcoblen. 

2* In both the forr.iulaticttis wages of handling equipment operators 
were not ccaasidered. In case of dynamic programming procedure 
ft'' 4#. 'hit- 'possible'; ^ti : 

■ prpgri0i^f:;«|jfroach 

\ V; ■fo'llifllJI 


If 
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W (J) 


Wago of opea^-ator of Jth type of handling 
equipaent 


then 


P 


% 


C (I,J,K) . X (I,J,K) + 2 ' W(J) ^Y (J,K) 


JtcP I«1 JfePI 

Seaae procedure can now bo used to so lye this problem. 


3, In case of dynanlo programing procedure, major fraction of the 
cocpitor tiao for finding solution at a grid point is taken in 
finding correi^ondlng block in the previous stage for state 
after trfsnsfornatim. Much of the computer time can be saved if 
an efficient method to search block for corresponding state can 
be devlsecl • 
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IDGIC PIDV/ DIAGRAM POR DHAHC 
DROGB/^MMING PROCEIORE 
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appmux - 5 , 

Data for tJia TOrioua problems ax'e given below 

problb: " 2 

(A) AGC®Ti 3 DE lOCiVflOIIS 

NEW rdOHIKE DOCAIION 1 2 


A 1 1 


Cb) ADTiMAlJB MillElEAlS MDIING systems 

um EXISIIIG H HI 

mCHBIE MACHINE 


(C) OTHM MTA 

. IDOMION mWIE,mG SYSTEM 
A f : . Q : . s ' s 


now Pi® A-?i 

Bqjiipaant 

1 

' . >18' 

•45 



Requirmmt 

2 

■■■ . ,21 ' 

,60 

- 


ilat«f:.lala ^'.j 

1 

2,227 

2,550 



haadling cost , 
la Rfci'llto'ia 

2 

■ " -3,275 ' 

5,750 


. . . *m 


1 

.18 

.51 


1.0 

S®<pir-«a«at 

2 



VA': 

1.0 

Materials ; 

1 




1640, 

haadillag cost 
la Ru/H&intt ! 

2 


2,880 


1640 

Purdaase oo»fe' 

1 



: V;'; 

24,000 

of 3aaadllag 



iliili 


24,000 

e^l|»eat -i 

2 
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jmuPiaxA*! .j H- 

rors ? 



tim 

i imm M Ki ii -A 1 * #»4 

I£!C/.‘i'I0H 

1 

■2 '3 , 

A 


1 

■t . -1 

B 


1 

0 1 : 

(b) kvmifysu lajTSim^ jOdling 

SYSIfflS 


Shf 

B 

¥1 

4if 

BII 

' 4k, «• 

P,Q,S 

- 

' P,Q,S 

B *«* ' . 

- 

P,Q 

P,Q,T 

(0) Ol’liER ' 

I£)G/*f ION 


: EANM 


A B 

■ B ' 

Q S T 

fhOl fi.TH B-K ■ 




Bq^ ipaiftt ■ 

1 0 

,18 


Ho^irasent . 

'■ 3 ; 0 ^ 

.21 

.60 - 

Materlfile Itaidling'' 
coet Is Rfc/toith 

1 0 

2,227 

2,550 • - 

■W 

3 0 






■'Ifnipaent ■ 

1 0 




3 0 

.18 



Material® Haadltog ' t 
Goat in Es*/llca3,tfe 




Purolmse cost 
of handling 
e^ipaent 




1,640 
1,640 
24,000 ,,, 
24,000 


■ •■■•I 





1 

2 

*15 

.45 

i.o 

Kuquiraut.nt 

1 

3 

.14 

.35 

3.0 


3 

1 

.14 

.35 

1.0 


3 

2 

.15 . 

.36 

1.0 


1 

2 

2,100 

2,400 

1,350 

Materials 

hundltof. 

1 

3 

1,890 

2,160 

900 

In 

3 

1 

1,890 

2,160 

900 


3 

2 

2, 100 

2,400 

1,170 


1 

2 

- 

- 

12,000 


1 




8,000 

BirchtsB© cost 

J ;. 




of hondltog 


1 


' *-«' 

8,000 

E{|ulpsB»nt 

■ 3 




, , 

•3 

' :/2 

■ - 


11,000 : 

fUOf A.-mi 

1 

' 2 ■ 

.24 

.54 

1.0 

. tpnwat 

1 

•■ 3 . 

' V '; . 33 ; 

.63 

1.0 

Requlroaent 

3 

1 

,21 

,48 

1.0 


3 


,24 

.54 

1,0 


1 


2,217 

3,223 

1,249 


1 

: , .5, ■ 

3,158 

4,012 

2,403 

.Materials ■ 



1,097 

' HaadUag ooat ; : 


1 

1,827 

2,730 

■ in Rs./teats : 

? ■ 


1,249 


i'wMi 

iil#i 

. 2,217 

3,223 


1 ■; 

iiiiii 



11,000 



Iflillils 



16,000 

]purciia®@ cost 
of bandliiag 

1 

'li i®i 

iilil 

i®if 



10,000 

©^Ipmint 

ilsi isis 

®Siil 


iiliSiiii 




Firot '3 j^low fCiihs are consida-ed as Problem - 3 . 


giOBlBM 5 end 6 

(a) Aocmism mciitiows 


KEW ».0EUIS 

ItOOATION 

1 

2 

3 

A 


1 

0 

1 

B 


1 

1 

1 

AETIEIIATI fMDMJO SYSTBIS 




A. 

B 

n 

FII 

Fin 

A 

*m ■ 

Q,T 

P,Q,S 


B P,Q 

' 0mk . 

B,Q,S 


P,Q,S 


(c) OTIffiR MA 


EOT WlTH -A-^m 
l^ipaeiat '1 


Material® 
handling coat 
in R4»/l6lm'tii 

'furcitoaae 

oo®t 


LOCATION 
A ' 3 


1 

5 

3 

1 


0 

0 

0 

0 

0 


HAIII][DIF& SYSTEM 

: : Q : ; S; 

.30 

.50 

800 

1,806 


60 


IIOI? 


Eqt ipDi4*^ i ■ 

1 

0 

♦40 

.70 

1.0 


Rj.'tp! h'ti .r,nt 

5 

0 

.40 

.70 

1.0 


Mtstcrinla 

1 

0 

841 

1,484 

563 

•rm 

hti’ndl kj/' asat 
in P,/t<!c;na; 

3 

0 

841 

1,484 

565 , 

- 

PurcitHjiij 

1 

0 


(N-. 

12,000 

, 

ouat 

3 

0 


. ,- ■ ■■ 

12,000 

- 

mm r/jH b-* 

1 

2 

.60 

.70 


■ 

E^iipaont 

1 , 

3 

.20 

.30 


- 

Sc^trsaant 

3 

r 

■ ,«:10' :, 

.10 


;- 

, I 

^'3 

■ z ' 

.40 

.50 


- 


1 

■ ■' 2 " , 

514 

761 



MBtWiol® 

1 

', ;3 ,' 

410 

475 



handltof ooBt 
to R«*/ltoto , 

3 

1 

171 

242 




[3.: : 



394 




■ mm ' 


MqxipRmt-:- 

Ee^l^WAl; 


K|‘f iiipf IIP- '^ ;• ji# V li: 'li' 


liiiiftiitiiMiite 




ilsSIlii 




M6t©ri®x» ;1 ,-i; •- 'Xi'-^r ■ '•>'■ '■ 
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1 
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- 

- 

12,000 


3 

2 

- 
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i 1. 

2 
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905 
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1 
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cost of 
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3 
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3 




o 

o 
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The first 5 Slow paths are considered as Problem - 5 
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Co) OTHER m.1% 


ABUmUBD HANiaiNG SYSTBf 
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/JPPMDIX « 4 

SAMPLE WTPlIf m DrmmG PROG®AMMIIf& HOOEUmE 

PROBOT - 2 

ifc sotipl© output for ProblcTi -■ 2, data for #iicla are giyea 
in Appendix « ”5 te Illustrated in thia /appendix. 


IMCRiME® IN milCflON - 1 » Rs, 1,000 

lICRaiBNf IN HHECflON - 2 » . .03 

IKOHMNT IN LCRECf LON - 3 « .03 


mNDLIlIG 

SrSTEM 


PLOW PAfH A-Fl 
I.OGI.TION C0P A - 1 


BI>OOK DlA,fl<mL POINTS 

i 2 , 


HANLUNG 

CGSa? 


(o, 0, 0) 

C§, 6, 15) 

99,999 

0 

(5, 0, 0) 

(60, 6, 54) 

'2,550 

2 

Co, 6, 0) 

". ,(5, 34,; 15) 


1 

Co, 0, IS) 

(5, 6, 34) 


IliiBliilil 

(60, 0, o) 

Cs, 6, o) 

Co, 6,;is) 

■•:;:;''C60,:;54fl|||| 
y-;(5,. .Ptiflill 

aiiiiisftiii 

..■■'■■.-a::'; ■ . 

2,227' 

A,?;,: w-J- ■i'-.'av -i-r;'. %-■■: ly-i'a-ia 
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Co, 0 , o) 

(5, 7, a). 

99,999 

0 

(S» 0, o) 

(60, 7, >1) 

3,150 

' 2 ■ 

Co, It o) 

(5, 54, 20) 

.3,273. '■ 

'I':' 

(o, 0, 2Q) 

C5,.7, 34) ■' 

3,750' 

' ■ 2' 

C&o, 0 , o) 

(85, 54 , 54) 

3,275 

: 1 '. 

ih 1$ 0) 

(60, 34, 34) 

3,275' 

- ' 1 

Co, 7, m) , 

.(5, 34, 34) 

3,275. 
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